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Abstract. We present new and improved evolutionary calculations for carbon-oxygen white dwarf (WD) stars 
appropriate for the study of massive ZZ Ceti stars. To this end, we follow the complete evolution of massive 
WD progenitors from the zero-age main sequence through the thermally pulsing and mass loss phases to the 
^ ' WD regime. Abundance changes are accounted for by means of a full coupling between nuclear evolution and 

time-dependent mixing due to diffusive overshoot, semiconvection and salt fingers. In addition, time-dependent 
^0 ' element diffusion for multicomponent gases has been considered during the WD stage. Emphasis is placed on the 

, chemistry variations along the whole evolution. In particular, we find that before the ZZ Ceti stage is reached, 

^ ■ element diffusion has strongly smoothed out the chemical profile to such a degree that the resulting internal 

' abundance distribution does not depend on the occurrence of overshoot episodes during the thermally pulsing 

, phase. The mass of the hydrogen envelope left at the ZZ Ceti domain amounts to Mh ~ 2.3 xlO"*^ M0. This is 

about half as large as for the case when element diffusion is neglected. The implications of our new models for the 
pulsational properties of massive ZZ Ceti stars are discussed. In this regard, we find that the occurrence of core 
overshooting during central helium burning leaves strong imprints on the theoretical period spectrum of massive 
Q ' ZZ Ceti stars. Finally, we present a simple new prescription for calculating the He/H profile which goes beyond 

^ , the trace element approximation. 

C/3 . 

' Key words, stars: evolution - stars: abundances - stars: AGB stars: interiors - stars:white dwarfs - stars: oscilla- 

tions 

X, 

1. Introduction mode oscillations in the DA WDs. ZZ Ceti stars are well 

■ ■ ■ ' known to pulsate in a narrow effective temperature (Tes) 

Variable DA white dwarf (WD) stars or ZZ Ceti stars be- ij^^erval ranging from 10700 K ^ T,s^ 12500 K. Numerous 
long currently to one of the best established and most important studies have been devoted to exploring the pul- 
extensively studied class of non-radial pulsatmg stars sational properties of ZZ Ceti stars, amongst them Tassoul 
(see Gautschy & Saio 1995, 1996 for a review). These (^ggQ)^ Brassard et al. (1991; 1992ab), Gautschy et 

hydrogen-rich, pulsating WDs exhibit multiperiodic lumi- (^ggg) Bradley (1996, 1998, 2001). 
nosity variations caused by gravity-modes of low harmonic 

degree (£ < 2) and periods in the range of 100-1200 s. Over the last past decade, the study of the pulsational 
Such modes were first thought to be excited by the k- pattern of variable WDs through asteroseismological tech- 
mechanism (Winget et al. 1981). However, the dominance niques has become a very powerful tool for probing the 
of convective energy transport in the driving region led internal structure and evolution of these stars. Indeed, 
BrickhiU (1991) and Goldreich & Wu (1999) to propose asteroseismological inferences have provided independent 
convective driving as the main driving mechanism for g- valuable constraints to fundamental properties such as 

core composition, outer layer chemical stratification and 
Send offprint requests to: L. G. Ahhaus stellar mass (Pfeiffer et al. 1996; and Bradley 1998, 2001 

* Member of the Carrera del Investigador Cientifico y amongst others). In particular, asteroseismology of mas- 
Tecnologico, CONICET, Argentina. sive ZZ Ceti stars has recently drawn the attention of re- 

** Fellow of CONICET, Argentina. searchers in view of the possibility it offers to place con- 

*** Fellow of CONICET, Argentina. straints on the crystallization process in the interior of 
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WDs (Montgomery & Winget 1999). This has been mo- 
tivated by the discovery of pulsations in the star BPM 
37093 (Kanaan ct al. 1992). a massive ZZ Ccti star which 
should be largely crystalhzed (Winget el al. 1997). 

Needless to say, a detailed WD modeling as well as a 
complete and self-consistent treatment of the evolution- 
ary stages prior to the WD formation are required for 
an adequate interpretation of the pulsational patterns. In 
this regard, the outer chemical stratification, the relics of 
the nucleosynthesis and mixing processes which occurred 
along the asymptotic giant branch (AGB) evolution (see 
D'Antona & Mazzitelli 1990), is a matter of the utmost 
importance as far as asteroseismology of ZZ Ceti stars 
is concerned. The construction of ZZ Ceti models based 
on full evolutionary calculations has recently been under- 
taken. Indeed, Althaus et al. (2002) have presented de- 
tailed models for a O.563-M0 ZZ Ceti remnant based on 
the complete evolution of an initially 3 Mq star, using 
models which take into account the chemical evolution 
during the WD regime caused by time-dependent element 
diffusion. The exploration of the pulsational properties of 
such models has been performed by Corsico et al. (2001; 
2002), who have found that mode trapping effects are con- 
siderably weakened as a result of the smoothness of the 
diffusion-modeled chemical profiles. These studies consti- 
tute an improvement as compared with most of the exist- 
ing research in pulsating DA WDs which invokes diffusive 
equilibrium in the trace element approximation to assess 
the shape of hydrogen-helium chemical interface. 

In this paper we extend the calculations presented in 
Althaus et al. (2002) to the case of massive intermediate- 
mass stars. As compared with Althaus ct al. (2002) major 
improvements in the treatment of the abundance changes 
have been made. In particular, we developed a time- 
dependent scheme that fully couples abundance changes 
due to nuclear burning, mixing processes and element dif- 
fusion. Time-dependent mixing due to semiconvection and 
salt finger is fully taken into account as well as exponen- 
tially decaying diffusive overshooting above and below any 
formally convcctivc region. In particular, the presence of 
some overshooting below the convective envelope diiring 
the thermal pulses has been shown by Herwig et al. (1997) 
to yield third dredge-up and carbon-rich AGB stars for 
relatively low initial mass progenitors (see also Ventura et 
al. 1999 and Mazzitelli et al. 1999). In addition, the oc- 
currence of extra mixing below the helium-flash convection 
zone during the thermally pulsing AGB phase is supported 
by recent evolutionary calculations (Herwig et al. 1999 
and Herwig 2000). These studies show that the resulting 
intershell abundances are in agreement with abundance 
determinations in hydrogen- deficient post- AGB remnants 
such as PC 1159 stars. 

The main aim of this work is to present the first results 
of detailed and complete evolutionary calculations appro- 
priate for the study of massive ZZ Ceti stars. To this end, 
we follow the evolution of initially 7.5- and G-M© stars 
from the zero-age main sequence through the thermally 
pulsing phase on the AGB to the WD regime. Attention is 



concentrated on the chemistry variations along the whole 
evolution. We emphasize in particular the role of time- 
dependent element diffusion on the chemical abundance 
distribution at the ZZ Ceti stage. The exploration of the 
pulsational properties of massive ZZ Ceti stars in the light 
of our new models is likewise within the scope of this in- 
vestigation. Next, in Sect. 2, we briefly describe the main 
physical inputs of the models, particularly our treatment 
for the chemical evolution and overshooting. In Sect. 3 we 
present the main evolutionary results. Pulsational results 
are discussed in Sect. 4. In that section we also present a 
simple new prescription for calculating the He/H profile, 
which constitutes an improvement over the trace element 
approximation. Finally, Sect. 5 is devoted to making some 
concluding remarks. 

2. Input of the models and evolutionary sequences 

Here we describe at some length the main characteristics 
of our evolutionary code. We restrict ourselves to the 
main updates in the macrophysics, particularly the 
treatment of the chemical abundance changes. 

General description of the code: The results presented 
in this work have been obtained with the stellar evolu- 
tion code LPCODE we employed in our previous works. 
The code has been developed at La Plata Observatory 
and it is described in Althaus et al. (2002) and references 
therein. For the purposes of the present paper, the code 
has been substantially modified particularly with regard 
to the treatment of the abundance changes, modifications 
which will be described later in this section. 

Briefly, the code is based on the method of Kippenhahn 
et al. (1967) for calculating stellar evolution. Envelope in- 
tegrations from photospheric starting values inward to a 
fitting outer mass fraction (close to the photosphere) are 
performed to specify the outer boundary conditions. The 
independent variable is ^= In (1 — Mr/M^) and the depen- 
dent variables are: radius (r), pressure (P), luminosity (Z) 
and temperature (T). The following change of variables is 
considered in LPCODE: 

p{n+l) ^ p{n) + In (1 + 

l{n+l) ^ (1) 

with ue, Up, Ux and ui being the quantities to be iterated 

AT AP Ar 

that are given by ue = Up = "r. = ^ 

and Ui = AZ, where superscripts n and n -I- 1 denote the 
beginning and end of time interval (see Kippenhahn et al. 
1967 for definitions). Here, ^ = In T, a; = In r and p = 
In P. Thus, the Henyey iteration scheme is applied to the 
differences in the luminosity, pressure, temperature and 
radius between the previous and the computed model. 
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As for the constitutive physics, LPCODE employs 
OPAL radiative opacities (including carbon- and oxygen- 
rich compositions) for arbitrary mctallicity from Iglesias 
& Rogers (1996) and from Alexander & Ferguson (1994) 
for the low-temperature regime. In particular, opacities 
for various metallicitics arc required during the WD 
cooling regime in view of the metallicity gradient that 
develops in the envelope of the models as a result of 
gravitational settling. The equation of state for the 
low-density regime comprises partial ionization for 
hydrogen and helium compositions, radiation pressure 
and ionic contributions. For the high-density regime, 
partially degenerate electrons and Coulomb interactions 
are also considered. For the WD regime, we include 
an updated version of the equation of state of Magni 
& Mazzitelh (1979). In LPCODE crystalhzation sets 
in when the ion coupling constant F = Z'^e^ / {r)kBT 
reaches the value 180. Neutrino emission rates and 
conductive opacities are taken from the works of Itoh 
and collaborators (see Althaus et al. 2002). A nuclear 
network of 34 thermonuclear reaction rates and 16 
isotopes has been considered to describe the hydrogen 
(proton-proton chain and CNO bi-cycle) and helium 
burning, and carbon ignition. Nuclear reaction rates 
are taken from Caughlan & Fowler (1988), except for 
the reactions i5N(p, 7)1^0, i5N(j3, a)i2C, i*^0(p, a)i-'^N, 
is0(p,7)19F, i2c(a, 7)160, i6o(a, 7)20Ne, i3C(a,n)i6o, 
i80(a,7)22Ne, ^"^neia^nf^Mg and 22Ne(a, 7)26Mg, 
which are taken from Angulo et al. (1999). In particular, 
the ^^G{a,^Y^O reaction rate given by Angulo et al. 
(1999) is about twice as large as that of Caughlan & 
Fowler (1988). 

To get a reasonable numerical accuracy, AGB models 
typically contained 1400 mesh points, except for the 
peak of the thermal pulses where about 2000 mesh 
points were required. Evolutionary time steps during the 
thermally pulsing phase ranged from a few days during 
the helium flashes and the subsequent phases where the 
third dredge-up may occur to some years during the 
stationary hydrogen-burning interpulse phase. Finally, 
mesh distribution is performed every three time steps. 

Chemical evolution: An important aspect of the present 
study is the modeling of the chemical abundance distribu- 
ti(m throughout all of the different evolutionary phases. 
To this end, we consider a time-dependent scheme for the 
simultaneous treatment of chemical changes caused by nu- 
clear burning and convective, semiconvective, salt finger 
and overshoot mixing. Needless to say, such a coupling 
between nuclear evolution and time-dependent mixing is 
much more physically sound than the instantaneous mix- 
ing approximation usually assumed in stellar modeling. 
In particular, this kind of treatment has been used by 
Mazzitelli et al. (1999) to study the lithium production 
by hot bottom burning in AGB stars (see also Ventura et 
al. 1999). 

In what follows, we present some details about the 
numerical method for the abundance changes included 



in LPCODE. Specifically, the abundance changes for all 
chemical elements are described by the set of equations 



if) ~ 



+ 



d 



dMr 



dY 



dMr 



(2) 



with Y being the vector containing the number fraction of 
all considered nuclear species. Here, mixing due to convec- 
tion, semiconvection, salt finger and overshoot is treated 
as a diffusion process which is described by the second 
term of Eq. (2) in terms of the mass coordinate Mr- The 
efficiency of mixing processes is described by appropriate 
diffusion coefficients D (see later in this section) . The first 
term of Eq. (2) gives the abundance changes due to ther- 
monuclear reactions, changes which, in the present work, 
are fully coupled to mixing processes. This term is ex- 
panded as a function of local abundances and cross sec- 
tions following the implicit scheme by Arnett & Truran 
(1969). After linearization, it can be written as 



Y 



At 



(3) 



where At is the time interval At = — and super- 
scripts n and n+l again denotes the beginning and end of 
time interval. Subscript j denotes the j^^ mesh point. F^ 
is a, N X N matrix {N being the number of nuclear species 
considered) and Aj a vector of TV dimension, both with 
elements involving abundances (Y^) and nuclear reaction 
rates at t". Using a three-point, finite-difference scheme, 
the second term of Eq. (2) can be approximated by 



d 



dMr 



dY 

dM^ 



n+l 



AY]+l+BY]+'+CY]+l 



(4) 



Here A, B and C are diagonal matrices with dimension 
N X N. Specifically, the non-vanishing elements of such 
matrices are given by 



(mj_i/2 - mj+i/2){mj-i - mj) 



{mj-1/2 - mj+i/2)(mj - mj+i) 



(5) 



(6) 
(7) 



where the subscripts j -I- 1/2 (j — 1/2) indicate an average 
for the mass shell between the mesh-points j and j + 1 (j — 
1). At the boundaries of the mass interval of integration, 
we have no mass flux {dY /dMr = 0), thus 

= (4-V)^I)|^. (8) 

bdry C'-'"- 



d 
'dMr 



{Anr^pfD 



dY 
'dMr 



A second-order Taylor expansion in the abundances 
leads us to 



^ {^^f{p\^D\,2 

{mi — 77x2)2 



an = 



(9) 
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at the top boundary, and to 



(m. 



7-1 







(10) 



at the bottom boundary (characterized by the J**^ mesh- 
point). 

Eqs. (2), (3) and (4) lead to the following system of 
linear equations to be solved simultaneously for the new 
chemical abundances at time 



where Tj = T^^ + - Bj and M.^ = Aj + V^/At. 
T^^ is a N X N diagonal matrix with elements^ 1/At. 
Fortunately, Eq. (11) has a special structure (tridiago- 
nal) which enables us to minimize storage of matrix coeffi- 
cients. Schematically, we have (Aj, J-'j and Cj are matrices 
and Y]+'^ and Mj are vectors) 



Ai C4 



\ 



Aj-3 Tj-2 Cj-i 

A J -2 J^j-i Cj 

Aj-i Tj ) 



^ 2 



rn+1 



' J-2 
' J-1 



( Ml \ 
M2 
M3 



Mj-2 

Mj-i 
\ Mj J 



(12) 



To solve this set of equations, we followed the method 
considered by Iben & MacDonald (1985). Eq. (2) is ap- 
plied to radiative {D = 0) zones and to convective, semi- 
convcctivc, overshoot and salt finger regions provided the 
diffusion coefficient D is appropriately specified for each 
process according to the adopted treatment of convection. 
We consider a total of 16 isotopes (TV = 16): ^H, ^H, ^He, 
^He, ^Li, ^Be, ^^C, ^^C, "N, i^N, isq, i^q, i^O, ^^F, ^ONe 
and ^^Ne. In LPCODE, abundance changes are performed 
after convergence of each stellar model (and not during it- 
erations). For a better integration precision in the nuclear 
evolution computation, each evolutionary time step is di- 
vided into 5 chemical time steps. 

^ Note that Eq. (11) couples nuclear evolution to the current 
composition change due to mixing processes. 



The evolution of the chemical abundance distribution 
caused by element diffusion during the whole WD evolu- 
tion constitutes an important point of the present work. 
In our treatment of time-dependent diffusion we have con- 
sidered gravitational settling and chemical and thermal 
diffusion for the following nuclear species: ^H, ''He, ''He, 
^^C, ^'*N and ^®0. The chemical evolution resulting from 
element diffusion is described, for a given isotope i, by the 
continuity equation as 



9t /diff 



d 



dMr 



(13) 



where Wi is the diffusion velocity. We have adopted the 
treatment of Burgers (1969) for multicomponent gases. In 
this work, our focus is on the chemical evolution occurring 
quite deep in the star, thus radiative Icvitation (which is 
expected to modify the surface composition of hot WDs) 
has been neglected. In terms of the gradient of ion densi- 
ties, diffusion velocities can be written in the form 



Wi = wf 



dlnrii 



ions(j) 



dr 



(14) 



where wf* stands for the velocity component due to 
gravitational settling and thermal diffusion. Details arc 
given in Althaus & Benvenuto (2000); see also Gautschy 
& Althaus (2002) for an application to DB WDs. Unhke 
our previous works, the abundance changes resulting from 
element diffusion are fully coupled to nuclear reactions. 
To this end, matrix coefficients in Eq. (12) have been 
appropriately modified. 

Overshooting: In the present study we have included time- 
dependent overshoot mixing during all pre- WD evolution- 
ary stages. The scheme for the abundance changes de- 
scribed above enables us a self-consistent treatment of 
diffusive overshooting in the presence of nuclear burn- 
ing. In particular, we have considered exponentially de- 
caying diffusive overshooting above and below any for- 
mally convective region, including the convective core 
(main sequence and central helium burning phases), the 
external convective envelope and the short-lived helium- 
flash convection zone which develops during the thermal 
pulses. Specifically, we have followed the formalism of 
Herwig (2000) based on the hydrodynamical simulations 
of Freytag et al. (1996), which show that turbulent ve- 
locities decay exponentially outside the convective bound- 
aries. The expression for the diffusion coefficient in over- 
shoot regions is 



Dos = Do exp 



-2 z 

"h7 



(15) 



where Dq is the diffusion coefficient at the boundary of 
the convection zone (see below), z is the radial distance 
from the edge of the convection zone, = fHp, where 
the free parameter / is a measure of the extent of the 
overshoot region, and Hp is the pressure scale height at 
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the convective boundary. In this study we have assumed 
/ = 0.015, which accounts for the observed width of the 
main sequence as well as the intershell abimdances of 
hydrogen-deficient post-AGB remnants (see Herwig et al. 
1997; 1999, Herwig 2000 for details; see also Mazzitelli et 
al. 1999). 

Treatment of convection: For this work, we included in 
LPCODE the extended mixing length theory of convec- 
tion for fluids with composition gradients developed by 
Grossman et al. (1993) in its local approximation as given 
by Grossman & Taam (1996). These authors have devel- 
oped the non-linear mixing length theory of double diffu- 
sive convection that applies in convective, semiconvective 
and salt finger instability regimes. According to this treat- 
ment, the diffusion coefficient D in Eq. (2) characterizing 
such mixing regimes is given by 

D = £a (16) 

where £ = aHp is the mixing length and a the turbulent 
velocity. The value of a is determined by simultaneously 
solving the equations for the turbulent velocity and 
flux conservation (Eqs. 9 and 17 of Grossman and Taam 
1996). In this work, the free parameter a is taken to be 1.5. 

Mass loss: Our treatment of mass loss is that of Blocker 
(1995). In particular, during the AGE evolution, the mass 
loss rate is given by 

M = 4.83 X IQ-^M^a'^s^'-'^R [MqIvt] (17) 

where Mr is the usual Reimer rate given by 
Mr = 4 X IQ-^^ri^LR/M with L, R and M the lu- 
minosity, radius and mass of the star in solar units. 
This formulation, based on dynamical computations for 
atmospheres of Mira-like stars, accounts for the strong 
increase of the mass loss rate and modulations by thermal 
pulses expected in the course of the AGB evolution. 
We set the value of the free parameter rjji somewhat 
arbitrarily at 1 in order to get a reasonable number of 
thermal pulses on the AGB. A more careful choice for 
the value of tjr is not relevant for the purposes of the 
present paper. In our computation, mass loss episodes 
taking place during the stages of core helium burning and 
red giant branch are considered according to the usual 
Reimer 's formulation with r?R=l. 

Global properties of the evolutionary computations: A 
primary target of this work is the construction of detailed 
massive DA WD models with cores appropriate 

for the pulsational studies. Because the self-consistent 
solution of nuclear evolution and time-dependent mixing 
demands a considerable increase of computing time, 
particularly during the AGB phase of the WD progenitor, 
we restrict ourselves to examining two cases of evolution 
for the progenitor: 



— As for the first case, which will be hereinafter referred 
to as sequence NOV, we computed the complete 
evolution of an initially 7.5-Mq stellar model without 
overshooting. The sequence was followed from the 
zero-age main sequence all the way from the stages of 
hydrogen and helium burning in the core up to the 
thermally pulsing phase at the tip of the AGB. After 
experiencing several thermal pulses and as a result 
of mass loss episodes, the progenitor departs from 
the AGB and evolves towards its WD configuration. 
Specifically, during the 7th pulse, an enhanced mass 
loss rate is invoked to terminate the AGB phase. 
Evolution was pursued to the domain of the ZZ Ceti 
stars on the WD cooling branch. The final mass of the 
WD remnant is about 0.936 M©. 

— As for the second case (sequence OV), we followed 

the evolution of a 6-M0 progenitor from the main se- 
quences until the star underwent four thermal pulses. 
Here, diffusive overshooting was applied during all evo- 
lutionary stages and to all convective boundaries. This 
sequence was done with the aim of assessing the in- 
fluence of overshooting on the WD internal chemical 
stratification and to explore its implications for the 
pulsational properties of massive ZZ Ceti stars as well. 
For this sequence, we have not computed the post- 
AGB evolution. Rather, we obtained an initial WD 
configuration by simply scaling the chemical profile 
built-up during the pre- WD evolution to the struc- 
ture of a hot WD model resulting from sequence NOV. 
Because we are interested in the final chemical strat- 
ification of ZZ Ceti models, this procedure is enough 
for our purposes. Here, the mass of the ^^C/^^0 core 
towards the end of thermal pulses amounts to 0.94 Mq. 

For both sequences, a solar-like initial composition 
(Y,Z)= (0.275,0.02) corresponding to Anders & Grevesse 
(1989) has been adopted. Interestingly, the mass of the 
resulting ^^C/^^0 core is quite similar in both sequences 
despite the very different initial mass values of the progen- 
itors (see also Mazzitelli et al. 1999). In fact, overshooting 
leads to larger convective cores during the core burning 
stages. We want to mention that the breathing pulse in- 
stability occurring towards the end of core helium burning 
has been suppressed (see Straniero et al. 2003 for a recent 
discussion about this point). Finally, for both sequences, 
time-dependent element diffusion was considered during 
the whole WD evolution. 

To the best of our knowledge, this is the first time that 
WD models appropriate for the study of pulsational prop- 
erties of massive ZZ Ceti stars are derived from detailed 
evolutionary calculations which include a self-consistent 
treatment of time-dependent element diffusion and nu- 
clear burning. We report below the predictions of our cal- 
culations, particularly for the chemical stratification and 
comment on their implications for the relevant pulsation 
properties. 
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— OV (overshooting) 

— NOV (no overshooting) 




I I I I \ I \ I \ I I 

4.4 4.2 4 3.8 3.6 3.4 

Log(V 

Fig. 1. The Hcrtzsprung-Russcll diagram for the evolu- 
tion of our 7.5- and 6-M0 sequences (NOV and OV, re- 
spectively) from the main sequence to the first thermal 
pulses. 

3. Evolutionary results 

3.1. Evolution of the white dwarf progenitor 

In this section, we describe the main results for the evolu- 
tionary stages prior to the WD formation. Attention will 
be restricted to analysing the relevant aspects for the WD 
formation, in particular the chemistry variations along the 
evolution. Although the full coupling between nuclear evo- 
lution and mixing implemented in LPCODE is appropri- 
ate for addressing problems such as hot bottom burning 
and lithium production in AGB stars, we will not explore 
them in this paper. Such specific issues would carry us 
too far afield, and we refer the K^ader to the works of 
D'Antona & MazziteUi (1996) and AlazziteUi et al. (1999) 
for details. Other recent relevant studies of AGB evolution 
are those of Sackmann & Boothroyd (1992), Vassiliadis 
& Wood (1993), Blocker (1995), Straniero et al. (1997), 
Wagenhuber & Groenewegen (1998) and Herwig (2000) 
amongst others. 

We begin by examining Fig. 1 which illustrates the 
Hertzsprung-Russell (HR) diagram for the WD progenitor 
from the main sequence to an advanced AGB phase. Solid 
and dashed lines correspond to sequences OV and NOV, 
respectively. A feature worthy of comment predicted by 
sequence OV is that central helium burning occurs mostly 
during a loop to the blue in the HR diagram. Indeed, the 
blue excursion begins when the central helium abundance 
by mass falls below « 0.75 and continues until the abun- 
dance has decreased below 0.1. For this sequence, the to- 
tal time spent during hydrogen and helium burning in the 
core amounts to 7.5 x 10'' yr, while for sequence NOV this 
time is 4.2 x lO'^ yr. Following the exhaustion of central 
helium and before the re-ignition of the hydrogen burning 
shell both of our sequences experience the second dredge- 
up. As a result, the surface composition is appreciably 




Fig. 2. Internal -^He, and ^'^O profiles for the 7.5-M0 
sequence (NOV) at three evolutionary stages after the end 
of core helium burning. A salt finger instability develops 
at Mr « 0.45 M0 (middle panel) which leads to a chemical 
redistribution there (bottom panel). See text for details. 



modified, particularly for sequence OV. In fact, for this 
sequence, the envelope hcliiim abundance rises from about 
0.30 (resulting from the first dredge-up) to 0.364. For 
sequence NOV, the second dredge-up brings the helium 
abundance to 0.328. It is also worth mentioning that dur- 
ing the second dredge-up phase, the mass of the hydrogen- 
exhausted core is strongly reduced. In particular, for se- 
quence NOV, the core mass is reduced from rs 1.50 Mq 
at the end of central helium burning to « 0.93 Mq after 
the dredge-up episode. For sequence OV, the hydrogen- 
exhausted core is reduced from 1.38 to 0.938 Mq during 
the dredge-up. 

The inner chemistry variations that take place along 
the red giant branch and early AGB evolution is well doc- 
umented in Figs. 2 and 3 for sequence NOV and OV, re- 
spcictively. Specifically, we show the evolution of the inter- 
nal helium, carbon and oxygen distribution as a function 
of mass for the evolutionary stages following the end of 
helium burning in the core. We begin by examining the 
results for sequence NOV depicted in Fig. 2. The upper 
panel shows the chemical profile when the helium con- 
vective core vanishes leaving a central oxygen abundance 
of 0.55 by mass. As evolution proceeds, the helium-rich 
layers overlying the former convective core are radiatively 
burnt, giving rise to an off-centred peak in the carbon and 
oxygen abundances. This is shown by the middle panel of 
Fig. 2, which corresponds to the moment when the star 
surface luminosity exceeds log(_L/LQ) = 3.72 for the first 
time after 4.25 x 10^ yr of evolution. Because the oxygen 
abundance, and therefore the mean molecular weight, de- 
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Fig. 3. Same as figure 2 but for the 6-M0 sequence (OV). 
Core oversliooting leads to a larger core by tlie 

end of central helium burning. 

creases inwards at Mr ~ 0.45 Mq, a salt finger instability 
characterized by a large diffusion coefficient develops at 
this point. The resulting salt finger mixing is responsible 
for the redistribution of the innermost ^^C/^®0 profile, re- 
distribution that takes place during the following 3 x 10^ 
yr and which is documented by the bottom panel of Fig. 
2. Note also that during this time interval the mass of the 
12(-i^16q (,Qj.g ]-^g^g incj-cascd considerably by virtue of he- 
lium shell burning, reaching 0.925 by the re-ignition of 
the hydrogen shell before the occurrence of the first helium 
thermal pulse. The behaviour for sequence OV is detailed 
in Fig. 3. Here, the size of the ^^C/^®0 core emerging from 
the convcctivc helium core burning (upper panel) becomes 
substantially larger (we remind the reader that sequence 
OV has a lower initial stellar mass than sequence NOV). 
Thus, a smaller fraction of helium remains unburned and 
so a shorter pre-AGB phase is expected. In fact, about 
4.2 X 10^ yr are needed to evolve from the end of core he- 
lium burning to the onset of the thermally pulsing AGB 
phase (bottom panel), as compared with the 9.3 x 10^ yr 
for sequence NOV. The carbon and oxygen distribution 
in the core is clearly different from the case without over- 
shooting, a feature which is expected to bear its signature 
in the theoretical period spectrum of ZZ Ceti stars (see 
Section 4). Note also that sequence OV is characterized 
by a somewhat larger central oxygen abundance (0.63 by 
mass) than expected for sciquence NOV. We want to stress 
again that breathing pulses have been suppressed in our 
calculations. The suppression of breathing pulses inhibits 
the formation of ^^C/^^0 cores with large central oxy- 
gen abundances (see Straniero et al. 2003 and references 
therein). However, when central helium has been substan- 
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Fig. 4. Upper panel: the time-dependence of surface lu- 
minosity (in solar units) for the 7.5-M0 sequence (NOV) 
during the first thermal pulses at the AGB. Bottom panel: 
the time-dependence of surface luminosity, hydrogen- (as 
given by CNO cycle reactions) and helium-burning lu- 
minosities (dashed and dotted lines respectively) for the 
same sequence as above. The time-scale is given in yr 
from the main sequence. A total of 7 thermal pulses have 
been computed before the WD progenitor departs from 
the AGB as a result of mass loss episodes. 

tially depicted, sequence OV experiences a small growth 
of the convcctivc helium core, which increases the central 
helium abundance from 0.07 to 0.10. This fact is, in part, 
responsible for the higher central oxygen abundance that 
characterizes sequence OV. 

Towards the end of the early AGB phase, hydrogen 
is re-ignited in a thin shell and the star begins to ther- 
mally pulse. Here, helium shell burning becomes unsta- 
ble (see Iben & Renzini 1983 for a detailed description of 
this phase). The time dependence of surface luminosity, 
hydrogen- and helium-burning luminosities for sequence 
NOV during the first five thermal pulses is given in Fig. 4. 
The interpulse period for this sequence is roughly 2.4 x 10^ 
yr. Note that the helium burning rate rises very steeply 
at each pulse. After experiencing 7 thermal pulses and 
considerable mass loss, the WD progenitor departs from 
the AGB and evolves towards the planetary nebula re- 
gion and eventually to the WD state. In our simulation, 
when departure from the AGB occurs, stationary helium- 
shell burning mostly contributes to the star luminosity, 
but shortly after, hydrogen burning takcis over. During 
the thermal pulses, mass loss episodes reduce the stellar 
mass from 6.15 to 0.936 M©. It is worth noting that during 
the interpulses the temperature at the base of the convec- 
tive envelope becomes high enough for hydrogen-burning 
reactions to occur. The nuclear processing at the base of 
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Fig. 5. Same as figure 4 but for the 6-M0 sequence (OV). 
Here, tlie first three thermal pulses are shown. 

the convective envelope is commonly referred to as "hot 
bottom burning", a process whose occurrence is strongly 
dependent on the treatment of convection. As we men- 
tioned, we will not discuss this aspect here (see D'Antona 
& MazziteUi 1996 and Mazzitelh et al. 1999 for recent 
studies); suffice it to say that for sequence NOV, temper- 
atures at the bottom of the convective envelope after the 
first five pulses reach about 6x 10^ K, which is high enough 
for an early, albeit moderate, onset of hot bottom burning 
to occur. In fact, the surface carbon abundance decreases 
from 0.0021 to 0.0019 during this phase. 

The presence of overshooting from below the convec- 
tive regions affects the evolution during the thermally 
pulsing AGB phase. The run of luminosities for sequence 
OV is detailed in Fig. -5. Note that the trend of the lumi- 
nosity evolution is different from sequence NOV. In partic- 
ular, the rise in the helium shell burning luminosity at the 
peak of the pulse is much more noticeable than in sequence 
NOV. This can be understood on the basis that overshoot 
from the bottom of the helium-flash convection zone car- 
ries some helium from the intershell region into deeper 
layers, where it is burnt at higher temperatures. This pic- 
ture can be visualized by examining Fig. 6, in which we 
show the internal abundance distribution versus mass co- 
ordinate during the second thermal pulse of sequence OV. 
The region illustrated comprises the base of the hydro- 
gen/helium envelope, the almost pure helium buffer, the 
intershell region and the top of the ^^C/^^0 core. The up- 
per panel corresponds to the moment at which the model 
is at the peak of its second thermal pulse. Here, the model 
is characterized by a short-lived (few years) convection 
zone, the product of the huge flux of energy caused by un- 
stable helium burning. This convection zone extends from 
Mr fa 0.9414 Mq down to the base of the helium-burning 
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Fig. 6. Chemical profiles for the 6-Mq sequence (OV) at 
three evolutionary stages during the 2'^'^ thermal pulse. 
At pulse peak (upper panel), the model is characterized 
by a "^He-rich buffer and an underlying region rich in ^He, 
^^C and ^^O (intershell). The ^He-flash convection zone 
and overshooting extending down to Mr ~ 0.938 Mq are 
responsible for the redistribution of the intershell abun- 
dances (middle panel). Note also that envelope convection 
has penetrated down to Mr= 0.9399 M0, sweeping away 
the "^He buffer and reaching ^^C-rich regions (third dredge- 
up). During the interpulse, the helium buffer is built-up 
again (bottom panel). Note also the formation of a small 
^^N-pocket at the base of the helium buffer. 

region at 0.9388 M0. At the bottom of this helium-flash 

convection zone there exists an overshoot region stretching 
down to the underlying carbon-rich layers (the ^^C/^^0 
core) at about 0.9383 Mq. As a result of this overshoot 
region, larger amounts of carbon and oxygen arc mixed 
up, as compared with the situation in which overshooting 
is neglected. This mixing episode is particularly impor- 
tant regarding the outer layer chemical stratification of 
the emerging WD remnant (see Herwig et al. 1997). As 
can be seen from the middle panel of Fig. 6, the result- 
ing abundances of the intershell zone are (^Hc/-'^^C/"'^^0)= 
(28/56/12) as compared with (74/23/0.4) for sequence 
NOV^. These results quantitatively agree with the predic- 
tions of Herwig (2000) for the case of lower stellar masses 
than considered here. Note also that envelope convection 
(and overshoot from below) has penetrated into deeper 
layers down to Af,.= 0.9399 Mq, reaching '^^C-rich regions 
(third dredgc-up'^). During this process, the helium buffer, 

^ The quoted ''He/'^^C values for the intershell abundances 
change somewhat with subsequent pulses, see Herwig (2000). 

^ A process which is also favoured by the occurrence of over- 
shoot at the bottom of the helium-flash convection zone. 
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which was built-up by hydrogen shell burning during the 
preceding pulse, is completely wiped out. These episodes 
take place while hydrogen shell burning is almost extinct. 
Sequence OV experiences important dredge-up of carbon 
even at the first pulse; however, hot bottom burning taking 
place during the following quiescent interpulse phase is so 
efficient (with temperatures at the base of the convective 
envelope exceeding 7 x lO'^ K) in converting carbon into 
nitrogen that the formation of a carbon star is avoided (at 
least during the first thermal pulses). Interestingly, about 
25 % of the luminosity of the star during this phase is 
produced within the convective envelope. 

Bottom panel of Fig. 6 shows that the helium buffer is 
built-up again during the interpulse phase and a next pulse 
is thus initiated. Finally, another feature predicted by our 
calculations is the presence of a small radiative ^^N-pocket 
at the base of the helium buffer (see also Herwig 2000 for 
a similar finding). In fact, during the third dredge-up dif- 
fusive overshoot has led to the formation of a small region 
in which hydrogen from the envelope and carbon from the 
intershell region coexist in appreciable abundances. When 
this region heats up enough to re-ignite hydrogen burn- 
ing, one of the main results is the formation of ^^N with 
abundances reaching about 0.5. However, the mass range 
over which the ^''N-rich region extends amounts only to 
PS 4 X 10~^ M0. Eventually, the ^^N-pocket is swept into 
the helium-flash convection zone during the next pulse. 

3.2. White dwarf evolution 

After considerable mass loss, the mass of the hydrogen en- 
velope is reduced so much that the WD progenitor aban- 
dons the thermally pulsing AGB phase and evolves to- 
wards large T^e values to becomes a WD. By the end 
of mass loss at Toff w 8000 K, the mass of the hydro- 
gen envelope of the 0.936-Mo remnant (sequence NOV) 
amounts to Mh « 5.5 xlO~^ M©. 200 yr later, the post- 
AGB remnant reaches the point of maximum Tdf and A'/h 
is reduced to 10^"'' Mq as a result of nuclear burning. We 
mention that a similar value for Mh has also been derived 
by Blocker (1995) for his most massive WD remnant. The 
time-dependence of Mh is displayed in Fig. 7. Notably, 
by the time the ZZ Ceti domain is reached, subsequent 
nuclear burning has further reduced the hydrogen mass 
to Mh « 2.3 xlO~^ M©. This occurs as the WD evolves 
through the T^a (logi/L©) range of 140000 : 25000 K 
(1.65 : —1.55). Most of this reduction is the result of nu- 
clear burning via CNO of hydrogen chemically diffusing 
inwards. For comparison, when diffusion is neglected, the 
final hydrogen content remains about 5.3 xlO~^ M© once 
nuclear burning via CNO becomes virtually extinct at T^s 
~ 50000 K after 1.6 x 10^ yr of evolution. We want to men- 
tion that because we have not invoked additional mass 
loss episodes during the planetary nebula stage or early 
during the WD cooling branch, the value of Mh should 
be considered as an upper limit. Because of the larger 
surface gravity characterizing massive WDs, the resulting 
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Fig. 7. Time-dependence of different luminosity contri- 
butions and mass of hydrogen envelope (Mh) (in solar 
units) for the O.936-M0 WD remnant of T.S-M© evolu- 
tion (sequence NOV): surface luminosity, Lgur, luminosity 
due to proton-proton reactions, Lpp, CNO bi-cycle, icNO: 
helium burning, Lhci neutrino losses, inou and rate for 
gravothermal (compressional plus thermal) energy release 
Lg. Thin lines depict Lqt^q, Lpp and Mh for the case when 
element diffusion is neglected. The approximate location 
of the ZZ Ceti instability strip is indicated with vertical 
dashed lines. Note that by the time the ZZ Ceti domain 
is reached, hydrogen burning is virtually extinguished. As 
a result of hydrogen chemically diffusing inwards, nuclear 
burning via CNO remains non-negligible at late times and 
the final hydrogen envelope mass becomes a factor two 
smaller than when diffusion is not considered. Time is in 
years counted from the end of mass loss. 

hydrogen envelope is less massive than in a typically 0.6 
MqWDs. For instance, for a 0.563 MqWD model, Althaus 
et al. (2002) find a hydrogen envelope mass of 7 x 10~^ 
Mq. 

In Fig. 7 we also show as a function of time the lumi- 
nosity contributions due to hydrogen burning via proton- 
proton reactions Lpp and CNO bi-cycle LcnOj helium 
burning Lne, neutrino losses Lnou, surface luminosity Lsur 
and gravothermal energy release Lg for the O.936-M0 
WD remnant (sequence NOV) from the end of mass loss 
episodes near the AGB to the domain of the ZZ Ceti stars. 
In Fig. 7 we also include the predictions for hydrogen burn- 
ing luminosities and hydrogen envelope mass for the situa- 
tion when element diffusion is neglected (thin lines). Some 
features of this figure deserve comment. In particular, ex- 
cept for first 200 yr when the remnant evolves to its WD 
configuration, nuclear burning never constitutes the main 
source of surface luminosity of the star. Note that had dif- 
fusion not been considered, nuclear burning via CNO cycle 
reactions would have ceased after only 10® yr of evolution. 
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a result which is in agreement with that of the Blocker 
(1995) calculations. By contrast, CNO reactions at the 
base of the hydrogen envelope remain non-negligible for 
a longer period of time (2 x 10* yr) in the WD evolution 
if diffusion is allowed to operate. This is because chem- 
ical diffusion carries some hydrogen inwards into the he- 
lium buffer and carbon upwards from the carbon-rich zone 
through the buffer layer (see later in this section), thus 
favouring the occurrence of nuclear reactions. However, by 
the time the ZZ Ceti domain is reached, hydrogen burning 
becomes virtually extinct. 

From the preceding section, we saw that the chemical 
stratification for both the core and outer layers of the WD 
remnant is quite different according to whether overshoot- 
ing is considered or not during the pre- WD evolution. We 
show the chemical abundance distribution at the start of 
the WD cooling branch in the upper panel of Figs. 8 and 
9 for the 0.936- and O.94-M0 WD remnants (sequence 
NOV and OV), respectively. The external chemical pro- 
file emerging from the thermally pulsing AGB phase (sec- 
tion 3.1) is characterized by a region rich both in helium 
and carbon (and oxygen if overshooting is considered), the 
relics of the short-lived mixing episode during the peak 
of the last helium pulse on the AGB, and an overlying 
helium-rich buffer. The mass of the helium buffer amounts 
to 3 X 10"^ and 9.4 x 10"'' M© for sequences NOV and 
OV, respectively. The buffer in sequence NOV is substan- 
tially less massive than in sequence OV simply because 
the WD progenitor in sequence NOV departs from the 
AGB at an early stage during the interpulse phase, whilst 
the initial chemical stratification for the WD in sequence 
OV corresponds to a pre- WD structure at an evolution- 
ary stage well advanced in the interpulse phase. Note also 
the signatures left by the third and last dredge-up episode 
during the last helium thermal pulse in sequence OV, in 
which convection penetrated into the carbon-rich region 
(see Fig. 6 for details). It is also clear from these figures 
the markedly different chemical profile resulting from the 
consid(Tation of core overshooting during central helium 
burning. 

On the cooling track, the abundance distribution of 
the WD remnant will be strongly modified by gravita- 
tional settling and chemical diffusion. This is illustrated 
in the bottom panel of Figs. 8 and 9 which show the chem- 
ical profiles at the ZZ Ceti stage for sequence NOV and 
OV. These figures emphasize the role of clement diffusion 
in the external chemical stratification of massive ZZ Ceti 
stars. Indeed, the shape of the chemical profile that was 
built-up during the AGB phase is virtually wiped out by 
diffusion processes acting during WD evolution. In partic- 
ular, near- discontinuities left by past mixing episodes in 
the shape of the external chemical interfaces (of primary 
importance in pulsation properties of WDs) are strongly 
smoothed out by the time the WD has approached the hot 
edge of the ZZ Ceti domain. It is worth noting that, ex- 
cept for the inner part of the core, element diffusion is so 
efficient that the resulting abundance distribution at the 
ZZ Ceti stage does not depend on whether overshooting 
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Fig. 8. Abundance by mass of ^H, "^He, ^^C and ^^O as a 
function of the outer mass fraction for the O.936-M0 WD 
remnant of T.S-Mq evolution (sequence NOV) at the start 
of the cooling branch (upper panel) and near the ZZ Ceti 
instability strip (bottom panel). Models are characterized 
by values (log(L/L0) , log Tcfr ) of (2.25, 5.27) and (-2.92, 
4.06) (upper and bottom panel, respectively). Clearly, el- 
ement diffusion strongly affects the chemical profile. 



during the thermally pulsing AGB phase is considered or 
not. In fact, the external chemical profile at the ZZ Ceti 
stage is quite similar for both sequences. By contrast, to- 
wards the innermost region, diffusion time scale becomes 
much longer than the evolutionary time scale, and the 
chemical profile remains therefore fixed during the whole 
WD evolution.^ Thus, overshoot episodes occurring dur- 
ing the pre- WD evolution leave recognizable features in 
the chemical profile of massive ZZ Ceti stars solely at the 
inner region of the core, features that are expected to leave 
their imprints on the theoretical period spectrum of these 
variable WDs (see section 4). As mentioned, our model 
is characterized by a chemical interface in which helium, 
carbon and oxygen in non-negligible abundances coexist, 
an interface which, at the ZZ Ceti stage, has extended 
appreciably as a result of chemical diffusion. 

The evolution of the chemical abundance distribution 
caused by diffusion processes can best be seen in Fig. 10, 
which illustrates the abundances by mass of ^H, ^He, ^He, 
i^C, "N and 1^0 for the O.936-M0 WD remnant (se- 
quence NOV) as a function of the outer mass fraction 
at various epochs characterized by values of logL/L0 and 
logTeff (numbers given in parentheses). Panel a shows the 



* Except for the minor mixing episode at the innermost 
i2q^i6q induced by an inward-decreasing molecular 



weight. 
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Fig. 10. Abundance by mass of ^H, ^He, ^He, ^^C, ^^N and ^^O in terms of the outer mass fraction for the O.936-M0 

WD remnant of 7.5-Mq evolution (sequence NOV) at four evolutionary stages characterized by values of logL/L© and 
logTefi (numbers given between brackets). Panel (a) corresponds to the start of the cooling branch, and panel (d) to 
the ZZ Ceti domain. At this stage, note the inward (outward) extent of the tail of the hydrogen (carbon) distribution 
resulting from chemical diffusion. 



chemical stratification at the start of the cooling branch at 
high luminosities. In the outermost layers, chemical com- 
position corresponds basically to that fixed by the second 
dredge-up episode during the early AGB phase. In the he- 
lium buffer, CNO abundances are different from those in 
the outer layers, because hydrogen burning in earlier evo- 
lutionary phases processed almost all the initial ^^C and 
^^O into ^*N. Rapidly, gravitational settling and chemical 
diffusion begin to play an increasing role. Indeed, gravi- 
tational settling causes hydrogen to float to the surface 
and heavier elements to sink down to such an extent that 
4.2 X 10'"' yr later and in the absence of any competing 
process, there are no metals in the outermost 10~^ Mq 
of the star (see panel b). With subsequent cooling, the 
effect of chemical diffusion is clearly seen at the chemi- 
cal interfaces where large abundance gradients exist. As 
a result, in the helium buffer there is a tail of hydro- 
gen reaching hotter layers and a tail of carbon diffusing 
outwards from the carbon-rich zone (panel c), thus caus- 
ing hydrogen burning via CNO cycle to be efficient for 
a longer period of time than when diffusion is neglected 



(see Fig. 7). Notably, the tail of the hydrogen distribution 
reaches a maximum depth by this epoch. Note also that 
by this time, diffusion has completely altered the shape 
of the chemical profile below the helium buffer. 7.3 x 10* 
yr later, the star reaches the ZZ Ceti domain (panel d). 
During this time, the tail of the carbon distribution has 
diffused outwards until reaching the outer boundary of 
the helium buffer. Also, the tail of hydrogen begins to re- 
treat outwards as a result of gravitational settling. This 
is so because electron degeneracy causes chemical diffu- 
sion to become less important and the inward diffusion of 
hydrogen is stopped. It is interesting to see that, despite 
the fact that the diffusion of hydrogen and carbon into 
the helium buffer zone has led to sufficiently large abun- 
dances of these two elements by the time the ZZ Ceti stage 
is reached, the temperature at the base of the hydrogen 
envelope is too low (w 5 x lO^K) for hydrogen burning to 
be of any importance. A last remark regarding the influ- 
ence of diffusion is related to the evolution of ^^C and ^^N 
abundances in the buffer layer. Indeed, initially the ^^N 
abundance in the helium buffer overwhelms that of ^^C 
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Fig. 9. Same as figure 8 but for the O.94-M0 WD rem- 
nant of 6-M0 evolution (sequence OV). Note that, except 

for the inner region of the core, clement diffusion acting 
during the WD evolution is so efficient that the resulting 
abundance distribution at the ZZ Ceti stage does not de- 
pend on whether overshooting during the thermally puls- 
ing AGB phase is considered or not. 



(see panel a), but before the ZZ Ceti stage is reached, 
is less abundant than ^^C throughout the star despite nu- 
clear burning having processed considerable ^^C into ^^N 
along WD evolution. The same finding, though to a lesser 
extent, has been reported by Iben & MacDonald (1986) 
for the case of a 0.6 MqWD. 

In what follows we explore the implications for the 
global pulsational properties of our models. 



4. Pulsational results 

For the pulsational analysis we have chosen two template 
ZZ Ceti models with Tcff w 12000 K, corresponding to 
sequences NOV and OV. In addition, we compare our 
results with the pulsational predictions for a O.94-M0, 
Teff = 12000 K WD model, as given by the stellar mod- 
eling considered in Montgomery & Winget (1999) (here- 
after MW model). We want to mention that despite the 
fact that the models are partially crystallized (between 
8% at the blue edge, to 25% at the red edge of 
the instability strip), we have computed the theoretical 
period spectrum assuming that the model interior is in 
a completely uncrystallized, fluid state. A complete dis- 
cussion of the pulsational properties taking crystallization 
self-consistently into account would be quite lengthy, so 
we have decided to postpone this to an upcoming com- 
munication, for which we will compute more massive WD 



models (and thus more suitable for the study of crystal- 
lized WDs) than attempted here. 

4.1. The Ledoux term and the Brunt-Vaisala frequency 

The Brunt-Vaisala frequency ( A'^) , a fundamental quantity 
in WD pulsations, is computed by employing the "modi- 
fied Ledoux" prescription. Specifically, A'' is given by (see 
Brassard et al. 1991 for details) 



^2 ^ 5^ Xt (y^^ - V + B) . 
P Xp 



(18) 



The Ledoux term B, for the case of a multicomponent 
plasma (M-component plasma), is given by 



B 
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EainXi 



i=l 



rflnP ■ 



(19) 



Here, Xt (Xp) denotes the partial logarithmic pressure 
derivative with respect to T (p) , V and Vad arc the actual 
and adiabatic temperature gradients, respectively, Xi is 
the abundance by mass of species i, and 



Xx, 



d\nP 



(20) 



P,T,{Xi^i} 



where P is the pressure. It is worth mentioning that the 
Ledoux term accounts explicitly for the contribution to N 
from any change in composition in the model (the zones 
of chemical transition). Brassard et al. (1992a) stress the 
relevance of a correct treatment of A'' at the interfaces of 
chemical composition transition zones in stratified WDs, 
particularly in connection with the phenomenon of mode 
trapping. 

In Fig. 11, the Ledoux term is plotted as a function 
of the outer mass fraction for NOV, OV and MW models 
(upper, middle and bottom panel, respectively). In ad- 
dition, we have plotted the abundance by mass of ^^O 
(dashed lines) and ^Hc (dotted lines) scaled by a factor 4 
for a clear and easy visual interpretation of the different 
features exhibited by B (see from Eq. 19 the strong depen- 
dence of B on the slopes of the chemical abundances). To 
begin with, note that the Ledoux term in models NOV and 
OV shows smoother peaks in the outer layers, as compared 
with the case of MW model. This different behaviour can 
be in part understood on the basis that MW invoke the 
diffusive equilibrium in the trace element approximation 
(see Tassoul et al. 1990) to calculate the chemical profile at 
the outer chemical interfaces, an approximation that leads 
to pronounced peaks in the B term. Note that the B term 
for OV model is characterized by a high-amplitude peak 
towards the innermost regions where the presence of over- 
shooting has led to a sharp variation of the ^^C/^^0 profile 
(see Fig. 3 and 9). This behaviour is in sharp contrast with 
the situation expected from the case in which core over- 
shooting is neglected (upper panel of Fig. 11). Note that 
the MW model also shows a peak in the core. However, 
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Fig. 11. The Ledoux term B in terms of the outer mass 
fraction as predicted by the NOV, OV and MW models 

(upper, middle and bottom panel). In order to make eas- 
ier the identification of the various features exhibited by 
B, the mass fraction of ^^O (dashed lines) and *He (dot- 
ted lines) is depicted, with their magnitudes arbitrarily 
increased by a factor of 4 for clarity. 

at variance with MW, models NOV and OV show an ex- 
tended bump at the bottom of the helium-rich zone, which 
again is a result of the broadness and smoothness of the 
diffusion-modeled profiles. 




Log(l-M/M,) 



Fig. 12. Same as Fig. 11, but for the logarithm of the 
squared Brunt- Vaisala frequency. 



In Fig. 12 we depict the logarithm of the squared 
Brunt- Vaisala frequency for the same cases shown in Fig. 
11. Each feature exhibited by the Ledoux term is clearly 
translated into N^. Except for the feature induced by the 
internal oxygen-carbon distribution, the Brunt- Vaisala 
frequency as predicted by the OV and NOV models is very 
smooth. In particular, note that the chemical transition 
regions lead to smooth bumps at logg —2.5 and -5.5. 
In contrast, the Brunt- Vaisala frequency corresponding to 
the MW model shows enhanced peaks at the location of 
all the chemical interfaces (see bottom panel of Fig. 12). 

4.2. Pulsational properties of the template models 

For the pulsation analysis of the NOV and OV models we 
have employed the same pulsational code as in Corsico et 
al. (2001; 2002). We refer the reader to those papers and 
references therein for details. The boundary conditions at 
the stellar centre and surface are those given by Osaki & 
Hansen (1973) (see Unno et al. 1989 for details). Pulsation 
computations for the MW model were carried out by em- 
ploying the same pulsational code as in Montgomery & 
Winget (1999)^. For each computed mode we obtain the 
cigcnpcriod Pk (being k the radial overtone of mode) and 
the dimensionless eigenfunctions 2/i , • • • , 2/4 (see Unno et 
al. 1989 for their definition). Following previous studies 
of WD pulsations, the normalization condition adopted is 
2/1 = 1 at the stellar surface. We also compute the os- 
cillation kinetic energy, (i^kin; see Eq. 1 of Corsico et al. 
2002) and the weight function, wf, given by Kawalcr et al. 
(1985). The weight function gives the relative contribution 
of the different regions in the star to the period formation 
(see for details Kawalcr et al. 1985 and Brassard et al. 
1992ab). Finally, for each model computed we derive the 
asymptotic spacing of periods as in Tassoul et al. (1990). 

For our template models we have computed adiabatic 
g- modes with £ = 1,2 and 3, with periods in the range 
expected for ZZ Ceti stars (50 s <^ Pk <, 1300 s). We 
begin by examining Fig. 13 which shows the values for 
the forward period spacing AP^ (= Pk+i — Pk) for £ = 1 
(upper panel) and i = 2 (lower panel) in terms of the 
periods computed. An inspection of the figure reveals that 
the amplitude of APk corresponding to the NOV model 
are typically lower as compared with the results of the 
MW model. This difTcrcncc is understood on the basis 
that the Brunt- Vaisala frequency of the NOV model is 
smoother than that of the MW model (compare upper 
and bottom panel of Fig. 12). The marked smoothness of 
the APk distribution in the NOV model implies that mode 
trapping is appreciably diminished in this model. 

In Fig. 14 we show the same APk-Pk diagram as in 
Fig. 13, but this time we compare the results of the OV 
model with the MW one. At first glance, the APk values 
for both set of computations exhibit similar amplitudes. 
However, a closer inspection of the plot reveals an inter- 

^ This code employ the outer boundary conditions as in Saio 
& Cox (1980). 
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Fig. 13. Period spacing values for I — 1 (upper panel) and 
2 (bottom panel) in terms of tlie computed periods, Pk for 
the NOV model and MW predictions. Filled triangles and 
solid lines (empty dots and dotted lines) correspond to the 
NOV (MW) model. 



esting feature. In fact, the APk distribution corresponding 
to the OV model clearly shows periodic minima with val- 
ues between minima almost constant and tending to the 
asymptotic value (not plotted). As we shall see towards 
the end of this section, the minima in APfc correspond 
to modes partially confined to the deepest regions of the 
OV model as a result of the pronounced A''^-peak showed 
by middle panel of Fig. 12. By contrast, the APfc distri- 
bution for the MW model shows clear signals of mode 
trapping with different amplitudes caused by the mode 
trapping due to the presence of all chemical interfaces in 
this model. 

Finally, in Fig. 15 we compare the APfc distribution 
corresponding to the NOV and OV models. We can clearly 
observe that the non-uniformity in the period spacing is 
much more apparent in the OV model, as expected from 
the shape of A^^ corresponding to these cases (Fig. 12). As 
stated before, the minima in the APfc distribution for the 
OV model are associated with modes partially confined to 
the high-density core region placed between the centre of 
the model and logg ~ —0.7. Consequently, these modes 
are energy-enhanced, so they must exhibit maxima in the 
£^kin distribution. To demonstrate this, we show in Fig. 16 
the oscillation kinetic energy distribution corresponding to 
the OV WD model, where we have labeled in particular 
the i = 2 modes with overtones k = 17, 19 and 21. Note 
from Fig. 15 that the £ = 2, = 19 cigcnmode corresponds 
to a minimum of APfc , but it has actually a E'kin value near 
a local maximum (Fig. 16). 
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Fig. 14. Same as Fig. 13, but for OV and MW models. 
Filled dots and solid lines (empty dots and dotted lines) 
correspond to the OV (MW) model. 
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Fig. 15. Same as Fig. 12, but for OV and NOV mod- 
els. Filled dots (triangles) correspond to the OV (NOV) 
model. 



To convince the reader that minima in APfc correspond 
to energy-enhanced modes, wc depict in Fig. 17 the ampli- 
tude of eigenfunctions yi and y2 (corresponding to radial 
and tangential displacements of matter, respectively) for 
i = 2, k = 17, 19 and 21 modes at the central region of the 
OV model. In addition, we show in the bottom panel the 
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Fig. 16. Kinetic energy values for ^ = 1 and 2 in terms of 
the computed periods, Pk for the OV model. 



weight function corresponding to these modes. An inspec- 
tion of this figure reveals that the amplitude of the eigen- 
functions as well as the weight function for the fc = 19 
mode are noticeably larger as compared with the neigh- 
boring ones {k = 17,21). Note also the evident distortion 
of the shape in y2 and wf for the modes at the loca- 
tion of the acute peak in the Ledoux term (logg ^ —0.7; 
middle panel of Fig. 10). This feature, which is a direct 
consequence of the abrupt fall of Xibq, acts to enlarge the 
magnitude of j/i, t/2 and wf of the fc = 19 mode, but di- 
minishes the amplitude of these functions in the case of 
modes with k = 17 and 21. 

4.3. Beyond the trace element approximation 

The trace element approximation has been a workhorse 
of chemical profile calculations for many years. Sadly, as 
these and other calculations have shown, it provides an 
inaccurate description of the transition zone and its as- 
sociated Ledoux term: it is too narrow, too peaked, and 
its first derivative is not continuous. In this section, we 
present a new prescription which retains the simplicity of 
the trace element approximation but which addresses its 
problems. 

From timescale arguments (e.g., Michaud & Fontaine 
1979) and numerical experiments, we expect the ''^He/H 
transition zone to be near diffusive equilibrium, so we are 
justified in using Eq. (A5) of Arcoragi & Fontaine (1980) 
for the equilibrium profiles. Traditionally, the trace ap- 
proximation results from solving this equation in the limit 
in which one of the two species is considered to be a trace. 
However, assuming complete ionization of both species, 
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Fig. 17. Eigenfunctions yi, y2 and the weight function wf 
(upper, middle and bottom panel, respectively) for modes 
with £ = 2 and k = 17 (dotted lines), fc = 19 (solid fines) 
and k = 21 (dashed lines). See text for details. 

we find that Eq. (A5) may still be integrated in closed 
form. For the '*He/H zone, we obtain 



-7*(1+7)*^(1 + 27)^^(1 + 67)" 

Pq 



(21) 



where 7 = nuc/nn is the ratio of number densities of 
the two species, P is the pressure, Pq is some reference 
pressure, and C is a constant. Additionally, if we assume 
a geometrically thin hydrogen layer (which is certainly 
valid), then the pressure can be directly related to the 
mass depth, so that this equation becomes 



C-^ =75(l+7)M(l + 27)-To(l-f 67)-^, 

where Qs = l-M^/Af*, qh = l-M^/M*, and C ; 
Since the mass fraction of He is related to 7 by 



^Hc ^^Hc 



47 



Ah ?iH + ^Hc nuc 1 + 47 



(22) 



1.764. 



(23) 



Eq. (22) gives us a smooth relation for the mass fraction 
of He valid from large abundances (ATho ~ 1, 7 3> 1) to 
low abundances (ATho <C 1, 7 ^ 1). 

Since the treatment given above offers an easy imple- 
mentation of diffusion profiles, we would like to see how 
it compares with self-consistently computed profiles. In 
Fig. 18 we compare the profile calculated using Eq. H22|l 
with that obtained from an evolutionary calculation as- 
suming self-consistent, time-dependent diffusion. We see 
that our parameterized profile provides a good represen- 
tation of the overall profile shape, and that it reproduces 
the maximum height of the Ledoux term to within ~ 20%. 
Since the numerical implementation of this prescription 
is only slightly more expensive computationally than that 
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Fig. 18. A comparison of the hydrogen mass fraction (up- 
per panel) and the Ledoux B term (lower panel) for the 
profile computed using Eq. |j22) (solid line), which assumes 
diffusive equilibrium, complete ionization, and an ideal gas 
equation of state, with that assuming full time-dependent 
diffusion (dotted line). 

based on the trace element approximation, we recommend 
Eq. (|22|l as a reasonable replacement in codes which use 
a pre-specifled functional form for the profile shapes. 

The main discrepancy between the two profiles shown 
in Fig. f 8 is most likely due to the neglect of electron de- 
generacy pressure which is implicit in the use of Eq. (A5) 
of Arcoragi & Fontaine (1980). Since the deeper C-O/He 
transition zone is in a more degenerate environment, as 
well as being much farther away from diffusive equilib- 
rium, a derivation along the lines which led to Eq. (|22|l is 
more difficult to justify; we therefore defer a discussion of 
a simple functional form for its shape until a future pa- 
per. Finally, we note that the above prescription for the 
"^He/H profile is aimed mainly at calculating the Brunt- 
Vaisala frequency. As such, it may need to be truncated 
at large depths in order to prevent unphysical situations 
such as a thermonuclear runaway due to H burning. 

5. Conclusions 

In this work we have computed new and improved evolu- 
tionary models for carbon-oxygen DA WD stars appropri- 
ate for the study of massive ZZ Ceti stars. In addition, the 
implications of our new models for the pulsational proper- 
ties of massive ZZ Ceti stars have been explored. To this 
end, we have followed the complete evolution of massive 
WD progenitors from the zero-age main sequence through 
the thermally pulsing and mass loss phases on the AGB to 
the WD regime. Attention has been focused on the model- 



ing of the chemical abundance distribution. In this regard, 
we developed a time-dependent scheme for the simulta- 
neous treatment of chemical changes caused by nuclear 
burning and mixing processes. Salt finger mixing, semi- 
convection and diffusive overshooting above and below 
any formally convective zone have been fully accounted for 
during the pre- WD regime. In this work, we have taken 
into account an extended mixing length theory for fluids 
with composition gradients. Also, time-dependent element 
diffusion for multicomponent gases has been considered 
during the WD evolution. 

An important aspect of the study has been to explore 
the implications of the occurrence of overshooting dur- 
ing the pre- WD evolution for the pulsational properties of 
massive ZZ Ceti stars. To this end, we restrict ourselves 
to examining two cases of evolution for the progenitor: se- 
quence NOV based on the evolution of a 7.5-M0 initial 
mass star in which overshooting was not considered and 
sequence OV based on the evolution of a G-M© star with 
overshooting. For both sequences, the mass of the result- 
ing core is quite similar (« 0.94 Mq). This has 
allowed us in principle to compare the pulsational prop- 
erties of carbon-oxygen massive ZZ Ceti stars having the 
same stellar mass but being the result of the evolution of 
progenitor stars with markedly different initial masses. 

As for the main evolutionary results for the WD pro- 
genitor, we mention: 

— Both of the sequences experience the second dredge-up 
episode during which the surface composition is appre- 
ciably modified from the assumed interstellar abun- 
dances. 

— A salt finger instability develops after the end of cen- 
tral helium burning, which is particularly noticeable 
in the case without core overshooting. This instability 
is responsible for the redistribution of the innermost 
12q^16q pj-Qfjig occurring before the end of the early 
AGB phase. 

— Sequence OV exhibits a sharp variation of the ^^C/^^0 
composition around Af,. « 0.72 M©, variation which 
leaves its signature on the theoretical period spectrum 
(see below). 

— During the thermally pulsing phase, the temperature 
at the base of the convective envelope becomes high 
enough for hot bottom burning to occur. This is par- 
ticularly noticeable for sequence OV. 

— The presence of overshooting below the bottom of the 
helium-flash convection zone gives rise to an enhance- 
ment of oxygen and carbon in the intershell region be- 
low the helium buffer, as compared with models with- 
out overshoot. This is in agreement with Herwig (2000) 
for the case of low- and intermediate- mass stars. 

— Sequence OV experiences the third dredge-up. At the 
end of this episode, a ^''N-pocket is built-up at the 
base of the helium buffer. 

For the WD evolution, our main conclusions are: 
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— The mass of tlie liydrogen envelope left at the ZZ Ceti 
domain amounts to Mh « 2.3 xlO~^ Mq. This is 
about half as large as for the case when element diffu- 
sion is neglected. 

— As a result of chemical diffusion, CNO reactions re- 
main non-negligiblc for a considerable period of time 
(2 X 10^ yr) during the WD evolution. However, by the 
time the ZZ Ceti domain is reached, hydrogen burning 
becomes completely unimportant. 

— Element diffusion strongly modifies the shape of the 
chemical profile during WD evolution. In particular, 
before the ZZ Ceti stage is reached, diffusion has 
already smoothed out the chemical profile to such a 
degree that the resulting external abundance distribu- 
tion docs not depend on the occurrence of overshoot 
episodes during the thermally pulsing phase on the 
AGB. 

Finally, the implications for the pulsational properties 
are: 

— Models with core overshooting (OV models) exhibit a 
much more pronounced non-uniformity in the spacing 
of consecutive periods, at variance with the situation 
without core overshooting. 

— OV models are characterized by the presence of modes 
which are "core trapped", i.e., modes with enhanced 
kinetic energy having eigenfunctions with relatively 
large amplitudes at the innermost, high-density regions 
of the models. 

— Pulsational properties are insensitive to the occurrence 
of overshoot episodes during the thermally pulsing 
AGB phase. This is due to the powerful action of ele- 
ment diffusion during WD evolution. 

— The period spectrum for our models is markedly dif- 
ferent from that of existing computations based on the 
use of the trace element approximation for assessing 
the shape of the chemical interfaces. 

— The trace element approximation need not be made 
in the solution of the equation of diffusive cquilibriiim 
given in Arcoragi & Fontaine (1980), which results in 
smooth profiles which are closer to the results of time- 
dependent diffusion. 

Wc conclude that the pulsational spectrum of massive 
ZZ Ceti stars as predicted by our new and improved evolu- 
tionary models turns out to be substantially different from 
that of previous research. In particular, because element 
diffusion strongly smoothes out the chemical profiles, the 
mode trapping caused by the outer chemical interfaces is 
notably diminished. As a result, the pulsational properties 
become very sensitive to the shape of the innermost chem- 
ical profile, i.e., the occurrence of core overshoot during 
the evolutionary stages prior to the WD formation. This 
could in principle be tested from observations of massive 
ZZ Ceti stars. We believe that this aspect would deserve 
an exploration also in the frame of less massive ZZ Ceti 
models than considered here. 



In closing, we judge that the evolutionary modeling 
presented in this work constitutes a physically sound and 
solid enough frame for exploring the pulsational properties 
of crystallized ZZ Ceti stars. We will address this aspect 
in a future communication. 
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